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Abstract 

Analytic solutions of the mean field evolution equations for an 
N-component scalar field with 0(N) symmetry are presented. These 
solutions correspond to rotations in isospin space. They represent gen- 
eralizations of the classical solutions obtained earlier by Anselm and 
Ryskin. As compared to classical solutions new effects arise because 
of the coupling between the average value of the field and quantum 
fluctuations. 
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1 Introduction 



A class of analytic solutions to the classical equations of motion of the non- 
linear sigma model have been constructed by Anselm and Ryskin[0] who 
conjectured the possible relevance of such solutions to high energy heavy ion 
collisions where a large number of pions are produced in the final states and 
a part of them may be described by a coherent state wave function. Subse- 
quently, Blaizot and Krzywicki||] presented another type of classical time- 
dependent solutions which possess a symmetry under the Lorentz-boost in 
the direction of colliding nuclei and used them to describe the soft pion emis- 
sion in ultra-relativistic nuclear collisions. Interest in these solutions grew up 
because of their possible relevance for the observation of a disoriented chiral 
condensate |^ which may be formed when the chiral symmetry, tem- 
porarily restored in hot hadronic matter produced by the nuclear collision, 
gets broken again by the rapid cooling of the matter by expansion. 

These classical oscillating solutions correspond, in the quantum language, 
to coherent states of quantum fields, rotating in the internal symmetry space 
of the fields. One may consider them as arising due to the spontaneous break- 
ing of the internal symmetry, similar to nuclear deformation and rotation. 
In the latter case, it is known that the nuclear rotation affects the nuclear 
deformation through the coupling of the collective motion of the deformed 
nuclear mean field to the motions of intrinsic nucleonic degrees of freedom in 
itp. It is the purpose of this note to show that similar phenomenon arises in 
the case of oscillating coherent quantum fields when their dynamics is treated 
in the mean field approximation. For this purpose we use the framework of 
the time-dependent mean field theory of the quantum fields which can be 
derived from the Gaussian Ansatz to the time-dependent variational wave 
functional of the functional Schrodinger equation]^, P]: 

= Af exp {i{n\(p — 1^)) 

X exp (^-{if - !f\{^ + iJ:)\(p - if)^ , 

where v'(x) = {x.\lp) is the coordinate of the quantum field and the G, S, 
(p, 7f define respectively the (time-dependent) real and imaginary part of the 
kernel of the Gaussian and its average position and momentum. Although 
this approach appears to sacrifice the covariance with respect to Lorentz 
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transformation, one can show that the resultant equations of motion can be 
rewritten in a covariant form and becomes equivalent to the ones obtained 
from an explicitly covariant formulation]^. 



2 Mean Field Equations 

For generality we consider an N-component scalar field (fai'^), a=l,. . . N with 
0(N) symmetry, characterized by a bare mass mp and a coupling constant A. 
The index a will be referred to in what follows as flavor or isospin index. In 
the mean field approximation the evolution of the mean value (f of the field 
(often called condensate) |^ 

is governed by the following equation of motion 



a + ml + ^Lf"^ + ^trS'(x, x) + x) 
6 



^(x) = 0. (1) 



where x= {t, x), i^{x) stands for the N-component vector ^^aix), and = i^l 
+ . . . + if'lf. Here it is implicit that the first four terms carry the NxN unit 
matrix. In the previous equation, 5" is a NxN matrix operator in flavor space 
which is related to the kernel of the Gaussian wave-functional 

S{x,x) = (x|G(t)|x) 

and the trace runs over the flavor indices. S{x, y) is the Feynman propagator 
which satisfies 

□^ + m^(x) S{x,y) = i6^{x - y), 
with the boundary condition S{x,y) = S{y,x)] it is symbolically written by 

S = ^r^ (2) 

□ + m^(a;) — te 

where the NxN mass matrix is 

m^(x) = ml + —Lp'^(x) + —tiS{x,x) 
6 6 

+ X ip{x) + -S{x, x). (3) 



3 



In this expression the symbol <f{x) x <f{x) denotes the NxN matrix whose 
(a, b) matrix element is (fa{x)ipb{x). The previous equations (|l]) - @ are non- 
linear because the motion of the condensate involves the mass matrix m^(x) 
which in turn involves the values of the condensate and of the propagator. 
A partial differential form of the mean field equations can be written down 



by introducing the so-called mode functions[10|. Note that the first three 



terms in equation (|TD correspond to the classical approximation considered 
by Anselm and Ryskin. The next two correspond to the contribution of 
quantum fluctuations whose effect is the object of the present study. 



3 Rotating Solutions 



Let us look for solutions of the previous equations by means of the following 
Ansatz 

(p{x) = U{x)ip^^^ = exp{i{q ■ x)Ty}Lp^^\ 

where = {u, q) and Ty IS cL generator of rotation in the subspace of flavor 
1 and 2: 
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Lf^^^ can be interpreted as the mean field in the rotating frame and we assume 
it to be independent of space and time, pointing in the direction of the first 
flavor: 

( Vo\ 

(0) _ 



V '■■ J 

We also introduce the propagator S'^^^ in the rotating frame by 

S{x,y) = U{x)S^''\x,y)U\y), 

with 

^(0) 



(4) 
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We assume that the NxN matrix M is time and position independent. This 
imphes that S^'^'> (x, x) has the same property. Indeed 



(2^ 



ip-{x-y) 



with 



{P^. + q^.ry) {p^' + q'^Ty) - + is' 

as can be checked by comparing the action of the two operators on a plane 
wave state exp(i/c • y). 

For the particular form we assumed for (p and S{x,y) the self consistent 
mass reads 



m^ix) = U{x) 



ml + ^ipl + ^tTS^'\x,x) 



+ ^(^(0) x(^(o) + ^5(o)(x,x) 



Using the formula 

U{x)du,U^{x) = 9^ - iq^Ty, 

it can be seen that the previous equation implies that our Ansatz actually 
solves the mean field equations provided the mass matrix M in the propagator 
satisfies 

= ml + ^^l + ^trS^'\x,x) 

+ ^(^(0) x(^(o) + ^5(o)(x,x). 

This NxN matrix equation defines the mass gap M. It will be referred 
to below as gap equation in the rotating frame. To have a closed set, the 
previous equations must be supplemented by the relation satisfied by the 
condensate (f^^^ in the rotating frame 

-q^ + ml + ^ipl + ^tTS^''\x,x) 



+^sr^{x,x) 



(fo = 0. 
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where we have use the assumption that the condensate in the rotating frame 
is position and time independent and takes the form (1). This equation has 
a trivial solution (fQ = having 0(N) symmetry when q = 0. The interesting 
solution however (at least when g = 0) is expected to have a non-vanishing 
value of the condensate. 



4 Gap Equation in the Rotating Frame 

Let us first show that the solution of the gap equation is a diagonal matrix 
which is furthermore a multiple of the unit matrix in the subspace a=3,. . . N 
i.e. = Ml 5 ah with Ma = for a=3,. . . N. 

Let us consider the propagator in the subspace of the first two flavors 
in the momentum representation. It can be written as the following 
2x2 matrix 

5(0) _ 1 fp^ + q^-M^ -2tp-q 



A V 2ip- q p^ + q^ - Mf , 
where A is the determinant 

A = + g2 _ ^2) ^ ^2 _ ^2) _ . ^)2_ 

In order to build the quantity S^^\x,x) we have to integrate S^^\p) over p. 
Since off-diagonal elements are odd functions of momentum, we see that the 
assumption of a diagonal mass matrix is compatible with the structure of the 
gap equation. 

The equations determining the four quantities characterizing the solution 
Mf, M|, /i^ and ipo are the following 

Ml = ml + l^l + ltTS^'){x,x) 

(5) 

+ l[Si'J{x,x)+6a,l^l\, 

while the condensate satisfies 

(m', - ^y^l - qA ^0 = 0. (6) 
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The diagonal elements of 5*'-°^ (x, x) have the following expression 




where for the first two flavors a=l,2 one has 



i{p^ + q^- Ml - M| + Ml) 



(7) 



(p2 ^q2_ Ml){p^ + g2 _ M|) - 4(p ■ q) 



2' 



where the masses appearing in this equation are supposed to carry a vanish- 
ingly small negative imaginary part. For the flavors a=3,. . . N one has 



a quantity which is manifestly flavor independent. 

These formulae complete the construction of the rotating solutions of the 
mean field equations for the sigma model. 

5 Structure of the Gap Equation 

The previous equations involve divergent integrals which need to be regu- 
larized e.g. by the introduction of a cutoff in momentum A. For g=0 one 
recovers the usual static mean field equations for the vacuum state. All cut- 
off dependence may, in principle, be removed by the proper renormalization 
of the mass and coupling constants. In the following analysis, we instead 
treat the 0{N) model as an effective theory retaining the cutoff dependence 
exphcitly. 

In examining the structure of the coupled integral equations (|) - (0) in 
the presence of non- vanishing g^, it is instructive to absorb into the masses 
by rewriting the effective masses of the first two flavors by A^^ = Mf — q^ 
and M.\ = M| — g^, while keeping other masses unchanged, M.\ = M^ for 
a > 3. Then we find, for the first two flavors. 



■I 



— jj? + ie'' 



Ml 



Ml 



ml-q' + ^^l + ^irS^'\x^x) + g-Sfh^.^) 
lb o 

ml-e + ^^l + \irS^'\x,x) + ^S^^{x,x) 
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where 





2' 



{p^ - Mi){p^ - Ml) - A{p ■ q) 



2 ' 



while the condensate equation becomes 



(• 



Ml - -ifl Uo - 0. 



One can see in this form that a part of the changes of the gap equation 
caused by non- vanishing q is to shift the bare mass parameter ml ml — q'^. 
We expect therefore that time-hke q^ causes an effect equivalent to increase 
\ml\, hence it leads to an increase of the amplitude of the condensate. This 
effect may be compared to that of the centrifugal force in ordinary rotation 
which would cause elongation of a non-rigid rotating body. Note that for 
time-like one can find a frame where the spatial component of q^ vanishes 
and the condensate rotates uniformly with angular frequency cu — \fq^. On 
the other hand, for space-like q^ there is a frame where the condensate be- 
comes static and oscillates spatially with the wavelength A = 27r/ In 
this case, the amplitude of the condensate tends to diminish from its value 
at g = 0. These effects are purely kinematical and arises also in the classical 
limit where the amplitude of the condensate is determined by 

- mo ^(/7o = (classical) . (8) 

The remaining effect of non-vanishing q is to introduce coupling between 
the first two fiavor states. This effect resembles that of the Coriolis force in 
nuclear rotation which introduces a coupling between single particle states 
with different orbital angular momenta in a deformed potential. This is a 
genuine quantum effect which does not show up in the classical approxima- 
tion. 

We now proceed to estimate the significance of these effects by perturba- 
tion theory for a small . 
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6 Perturbative Calculation of g^^-dependence 

Let us denote by the change in when going from a square momentum 
zero to q^. The corresponding change in the diagonal matrix element of S^^^ 
may be expressed as 



5S^^}{x,x) = JJMl + K,q^ 



(9) 



The coefficients Ja and Ka can be obtained by using the series expansion of 
the propagator (0) and integrating out over the momentum p. The result is 
given by 



K3 



G\Ml), J^ = J, = ... = J^ 
■ ■ ■ = iTjv = 0. 



where we have introduced the following integrals 

d^p i 



G{m^) 
ml Ml) 



-Ml 



(27r)^ — m? + is 
d^p 



(27r)Mp2-M2)2(p2_M2)- 



and defined the function G'{im?) as the derivative of G{m'^) with respect to 
m^. Here Mq stands for the value of Mi when g = while the unperturbed 
masses of all other flavors [a = 2, . . . , N) are set to fi. They satisfy the gap 
equations at g = 0. 

The integral G{'m'^) is divergent and requires regularization. Using a 3- 
dimensional regularization in momentum space with a cutoff A it is found to 
be 



G{m^ 







" 

m 


- log 


(-) 


+ - 




\m J 





Differentiating with respect to the square mass gives 
1 



log (— 

\ m 
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On the other hand, the integral I{M^, M^) is finite. Using the standard 
integral representation of propagators 0], the following closed expression 
for this integral can be obtained 



Mi 



167r2 M2 - M2 



X 



1 



Mi 



Ml - Mi 



The linearization of the gap equation leads to the following relations for 
the mass changes in the first two flavors 



5 Ml = ^5^1 + ^tr55(°)(x,x) + ^ \ji5Ml + K^q' 



and 



6M^ 



3 L 



For other flavors one flnds 
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A closed system of linear equations can be constructed by including the 
linearized form of the evolution equation for the condensate (§) 



3Ml - -5^1 - q' 



0, 



(10) 



and the expression for the total change in the trace of the propagator 
tr(5^(°) (x, x) = Ji5Ml + Ja [5 Ml + {N - 2)5 n^' 

+ {Ki + K^)q\ 

The solution of this linear set gives for the changes in the masses 



and 



5Ml = (C162 - &1C2) 



5fi^ = (aiC2 — ciaa)- 
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The constants appearing in these equations are 
ai = -Jl, 02 = 1 + -Jl, 
h = ^-^J2, 62 = (iV+ 1)^^2-2, 



Cl = 1 - -Ki, C2 = 1 - -Ki - - — . 

3 3 3 Oi 

The quantity A in the above equations is the determinant of the hnear system 



/\\2 ^ 
A = (iV + 2)(-j J1J2-AJ1 + -J2-I. 



The change in the condensate is given by equation (R_ 



7 Discussion 

In the above resuh all quantum effects are contained in the parameters Ja 
and Ka which represent the changes in the vacuum fluctuation caused by 
non-vanishing q. The classical relation (|^) can be obtained from this result 
by setting Ja = Ka = 0. In this limit one finds 

—5(fl = 2q^ (classical) . (11) 
3 

To obtain an estimate of the quantum fluctuations we have chosen Mq to 
be of the order of the sigma mass i.e. about 500 MeV and set = 4. 
Assuming the value of the condensate Lpo to be the pion decay constant 
/tt = 93 MeV, the condensate equation at g = gives a coupling constant 
A = 86.7. For a momentum cutoff A = 1 GeV, the gapequation gives /i = 224 
MeV. Note that /i is a variational parameter characterizing the significance 
of the quantum fluctuation of pion fields; it should not be taken as the 
pion mass. For the previous values of the masses one has Ji=-. 00489, J2=- 
.0141, i^'i=. 000948 and K2= 0.00800. Ahhough the coupling constant is 
large compared to unity the constants Ki and K2 are sufficiently small to 
have Cl and C2 very close to unity. Actual values are 

Cl = .973, C2 = .812, 
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which means that the most important explicit dependence is the one arising 
in the equation for the condensate. The above values give the following 
expressions for the variations of the masses with momentum 



5Ml = %.l2q^ and = 1.48g^ 



The change in the condensate is found to satisfy 



-5^1 = 7.15gl 



(12) 



An increase of the angular velocity (time-like g^) thus produces an increase 
of the chiral radius, both in the classical approximation and in the mean-field 
picture, in agreement with the qualitative argument given above. The change 
is more rapid in the quantum case. Similarly for a space-like a vanishing 
value of the condensate is reached more rapidly in the quantum case. This 
is expected because quantum fluctuations smear out the effective potential 
and make symmetry breaking more difficult to reach. In the quantum case 
a transition occurs when 



At this point the mass of the first flavor vanishes which implies that a solution 
with a non zero condensate can no longer be obtained beyond this point. 
The corresponding classical excitation energy density of the meson field is 



It is interesting to compare this critical momentum for space-like conden- 
sate to the corresponding value in the classical limit obtained from (P) for 
(po = 0: 



where we have used the relation = —2rn^ in the tree (classical) approx- 
imation. The coupling to the quantum fluctuation thus works to suppress 
the appearance of static condensate with longer wavelengths. This implies 
that the appearance of the static chiral condensate is more sensitive to the 
size of the system than the classical treatment predicts. This would have 
important implications for the dynamics of chiral condensate in high energy 
nuclear collisions. 



= ql ^ Ml/l. 



(13) 



ql fl/2 ^ (123 MeV)^ 



ql = Ml/2 



(classical) 
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It is worthwhile noting that the analogy with nuclear deformation is lim- 
ited by the fact that in the nuclear case deformation disappears in the clas- 
sical limit contrary to symmetry breaking in the 0{N) model. Evolutions 
with angular velocity are however similar. 

In conclusion we have shown that, as compared to analyses based on clas- 
sical equations, new effects arise in a quantum framework as a result of the 
coupling between the evolutions of mean values and fluctuations. Our re- 
sult indicates that the coupling to the quantum fluctuation suppresses static 
condensate with spatial oscillation while the condensate oscillating in time 
is ampliflcd by the quantum effects. 
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